Cooling of nanomechanical resonator by thermally activated single-electron transport 
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We show that the vibrations of a nanomechanical resonator can be cooled to near its quantum 
ground state by tunnelling injection of electrons from an STM tip. The interplay between two 
mechanisms for coupling the electronic and mechanical degrees of freedom results in a bias-voltage 
dependent difference between the probability amplitudes for vibron emission and absorption during 
tunneling. For a bias voltage just below the Coulomb blockade threshold we find that absorption 
dominates, which leads to cooling corresponding to an average vibron population of the fundamental 
bending mode of 0.2. 
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Remarkable steps are now being taken towards achiev- 
ing conditions under which quantum effects are exper- 
imentally accessible in nano-electromechanical systems 
(NEMS)[l[. This is encouraging both for the prospects 
of realizing a plethora of applications that depend on 
our ability to control and monitor the coherent dynamics 
of nanometer-scale mechanical systems and for shedding 
light on purely fundamental issues, such as the nature of 
the crossover from classical to quantum physics 0. 

A number of proposals have been put forward in or- 
der to reach the extremely low temperatures T of order 
huj/ks, where quantum effects become observable in me- 
chanical resonators of eigenfrequency w. Their aim has 
been to replace the conventional dilution refrigerators by 
more efficient active cooling methods. 

Most of them are based on well-established principles 
for laser cooling of atoms and molecules Q, but alterna- 
tive approaches have been proposed by several authors 
who suggest that the coupling between mechanical and 
electronic degrees of freedom can be exploited for induc- 
ing energy to flow from the former to the latter All 
these schemes for ground state cooling are based on quan- 
tum resolved sideband transitions between discrete quan- 
tum levels of the refrigerant, a cooling mechanism that re- 
lies on energy conservation in order to suppress processes 
which involve emission of vibrational energy quanta (vi- 
brons) with respect to those that lead to absorption of 
such quanta. 

In this Letter we suggest a new mechanism for ground 
state cooling of a nanomechanical resonator, which is 
based on passing a current through the resonator under 
conditions such that the probability amplitude for tunnel- 
ing electrons to emit vibrons is much lower than it is to 
absorb them. As such it has the advantage that it does 
not require the refrigerant to have a discrete energy spec- 
trum, which puts fewer constraints on the experimental 
design. 

To be specific we consider the model system sketched 
in Fig. [I] where electrons are injected from the tip of a 




Figure 1. Sketch of the model system considered. A metal- 
lic carbon nanotube is suspended over a trench between two 
grounded metallic leads, while an STM tip placed a distance 
h above the nanotube is biased at a negative voltage —V. 



scanning tunneling microscope (STM) into a suspended 
metallic carbon nanotube. Low-temperature tunneling 
spectroscopy studies on a similar device have shown 
that inelastic electron tunneling can create a non-thermal 
equilibrium population of vibronic states in the nan- 
otube. Below, we will show that the probability for 
vibron emission can be suppressed as a result of the 
interference between two different mechanisms for cou- 
pling the mechanical and electronic degrees of freedom 
of the system. One of these mechanisms is the nanotube- 
position dependent probability amplitude for electron 
tunneling from the STM tip to the nanotube, the other is 
the electrostatic force on the nanotube when it is charged. 

It turns out that the effect of the interference depends 
on the voltage bias between the STM tip and the leads. 
Our analysis shows that the destructive interference is 
maximal for a bias voltage slightly below the threshold 
voltage for lifting the Coulomb blockade of electron tun- 
neling through the system. If the nanotube is weakly 
enough coupled to the environment, the suppression of 
vibron emission is strong enough to drive the nanotube 
to near its vibrational ground state and hence effectively 
"cool" the mechanical degrees of freedom. 

In order to analyze the dynamics of the nanotube and 
of the tunneling electrons in the quantum regime we in- 
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troduce a model Hamiltonian, 



H — He + H,n + Ht + He 



where 
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+ h.c, (4) 



and where a^lo and c"q' are annihilation (creation) oper- 
ators for electrons in the STM tip [a ~ S), in the leads 
(a = L), and in the nanotube, respectively. 

The first term of (P), He, describes the STM tip, the 
leads, and the nanotube as reservoirs of non-interacting 
electrons. The second term, Hm, describes the nan- 
otube's mechanical degrees of freedom, which we restrict 
to the fundamental bending mode considered as a simple 
harmonic oscillator with angular frequency w, being 
the annihilation (creation) operator for an elementary 
excitation (vibron) of this mode. 

Electron tunneling between the STM tip and the nan- 
otube and between the nanotube and the leads are de- 
scribed by Ht, the third term of the Hamiltonian, in 
terms of the tunneling amplitudes ts and t^. Here the 
operator e**^ changes the number N of excess electrons 
on the nanotube by one, e~^'^Ne^'^ = N + 1. Since ^5 
depends on the overlap between electronic states in the 
STM tip and the nanotube, it depends on the deflection 
of the tube through the operator X — Axgs{b'' +b), where 
Axgs = y^h/{2Mijj) is the displacement uncertainty in 
the vibrational ground state and M is an effective oscil- 
lator mass. For simplicity we assume that the STM tip 
is positioned above the midpoint of the nanotube (see 
Fig. [1]) and model the deflection dependence of the tun- 
neling amplitude as ts{X) = ts exp(X/A), where A is the 
characteristic tunneling length of the barrier (A ~ 10~^° 
m). This dependence amounts to a coupling between 
the electronic and mechanical degrees of freedom that 
we will refer to as a tunneling electromechanical (TEM) 
coupling. In contrast, the distance between the nanotube 
and the leads is fixed, so that tL does not depend on the 
nanotube deflection. 

The last term. He, in ([TJ describes the electrostatic 
interactions in the system, which we will treat in the 
framework of the capacitance model. In this approxima- 
tion He only depends on the total charge on the nan- 
otube and on the voltages applied to the bulk electrodes. 
Assuming the supporting leads to be grounded and that 
a negative electrostatic potential —V {V > 0) is applied 
to the STM electrode, we restrict our analysis to the 
Coulomb blockade regime in which at most one extra 
electron may reside on the nanotube. Under such condi- 



tions He can be written as [ 
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where Cs = Cs + Cg, Cs is the mutual capacitance be- 
tween nanotube and STM tip, Cg is the total capacitance 
between nanotube and ground, Vc = s/'2Cg is the thresh- 
old value of V for lifting the Coulomb blockade and — e 
is the electronic charge. 

In general, Cs and Cg both depend on the geometry 
of the system and therefore on the nanotube deflection. 
Here we will only take the dominant deflection depen- 
dance of the STM-nanotube capacitance into account. 
Hence Cs = Cs{h — X), where h is the distance be- 
tween the STM and the straight nanotube. For small 
displacements of the nanotube we may linearize the in- 
teraction Hamiltonian ([5]) and use an approximation that 
for Cs Cg takes the form 



He = Uc(V)N - SXN - a{X)V^ 



(6) 



where S = 2{dCs/dx)oVc5V , and 5V = Vc ~ V. The 
first term of ^ determines the Coulomb blockade effect 
in the absence of nanotube deflections, while the sec- 
ond is a deflection-dependent electromechanical interac- 
tion term. Due to a formal analogy with the interaction 
term in the model Hamiltonian for the polaron problem, 
we will refer to the origin of this term as a polaronic 
electromechanical (PEM) coupling. The last term of (O 
is a contribution that does not depend on whether the 
nanotube is charged or not. 

It is important for what follows that the sign of the po- 
laronic force constant 5^ in ([S]) depends on the bias volt- 
age. If the bias voltage is below the Coulomb blockade 
threshold, so that only thermally activated transport is 
possible, i.e. if SV > 0, then > and hence if charged 
by an electron the nanotube will be attracted to the STM 
tip. On the other hand, if SV < 0, then < and the 
charged nanotube is repelled from the STM. 

Note that the possibility to change the direction of the 
force by varying the bias voltage crucially relies on the 
discrete nature of the tunneling charge. If this charge 
could be arbitrarily small, then Vc — >■ and hence 

cx —V. For any (positive) value of V the polaronic 
force would therefore be negative and push the charged 
nanotube away from the STM tip, decreasing the tunnel- 
ing matrix element 0. 

As we have seen above, the electromechanical inter- 
action is described by two separate terms in the Hamil- 
tonian, one due to what we call TEM coupling and the 
other due to PEM coupling. The cooling mechanism to 
be discussed below results from the interplay between 
these different types of coupling. In order to analyze this 
interplay, it is convenient to apply a unitary transforma- 
tion that removes the polaronic term from the Hamilto- 
nian and instead makes the tunneling amplitudes depen- 
dent on both the midpoint position X of the nanotube 
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and its conjugate momentum P = ih{b^ -'b)/2Axgs. This 
is achieved by the transformation H H — UHU\ 
where U = eyip{iAxePN/h). Here Axe = S/2Mlo'^ 
is the difference between the equihbrium positions of 
the charged and neutral nanotube. To leading order in 
the small dimensionless parameters et = Axgs/X and 
£p = Axe/Axgs {Sp ~ 0.1-0.01) the transformed tun- 
neling Hamiltonian (U) is 

Ht = ts^{l - {et + ep)b + (e* - ep)b'') cjafc,s 

k,q 

+ tL'^{l - Spb + Epb'') alj^Cq + h.c. (7) 

k,q 

From Eq. (jT)), it follows that in the Born approximation 
the rate of inelastic single-electron tunneling from the 
STM tip to the nanotube accompanied by the absorption 
or emission (— ) of a vibron is 



rs,± = r5(e?+£2±2et£p) 



(8) 



where Ts = Ts{V,T) is the rate of elastic electron tun- 
neling across the STM-nanotube junction. The first (sec- 
ond) term of (jS]) gives the probability for tunneling as- 
sisted by either absorption or emission of a vibron due to 
the TEM (PEM) coupling alone, while the third term cor- 
responds to the "interference" between these two mecha- 
nisms in the case of vibron emission (— ) and absorption 
(-I-). Clearly, the probability for vibron-assisted electron 
tunneling is different depending on whether a vibron is 
absorbed or emitted and the difference can be controlled 
by the bias voltage since Sp oc Axe{SV). 

In particular, Axe > if SV > so that the inter- 
ference is destructive (constructive) for tunneling accom- 
panied by vibron emission (absorption). If SV < the 
situation is reversed in the sense that Axe < and the 
interference is constructive (destructive) for emission (ab- 
sorption) processes. 

The case of constructive interference for emission pro- 
cesses has been analyzed in Ref. 8, where it was shown 
that a promotion of emission over absorption processes 
may lead to an electromechanical instability of the sys- 
tem if V exceeds a certain dissipation-dependent thresh- 
old. Here we will focus on the reverse situation. 

A complete suppression of the emission processes 
would eventually drive the mechanical subsystem to its 
ground state. However, two more types of electronic tran- 
sitions that may generate vibron emission remain to be 
considered. The first is the tunneling of an electron from 
the nanotube to the STM. By virtue of time reversal 
symmetry, the mechanism responsible for the suppres- 
sion of vibron emission during tunneling from the STM 
to the nanotube stimulates the emission of vibrons dur- 
ing tunneling in the reverse direction. In order to make 
the effect of such transitions negligible in the energy bal- 
ance for the mechanical subsystem, an electron that has 



tunneled from the STM should escape from the nanotube 
to the leads before it can tunnel back to the STM by an 
inelastic transition. This requires that \ts\ ^ \tL\ and 
ksT <C eVc, where the latter constraint ensures an ex- 
ponential suppression of the probability for electrons to 
tunnel from the leads to the nanotube. These conditions 
are satisfied in the real experimental situation. 

In addition to the "backward" tunneling transitions, 
vibrons can also be emitted when electrons tunnel from 
the nanotube to the leads, but then only by virtue of the 
polaronic coupling mechanism (see Eq. (O). 

From Eqs. ([8]) and (O and the definitions of St and £p, 
it follows that the ratio beween the total rate of vibron 
emission and the total rate of vibron absorption reaches 
an absolute minimum for the bias voltage V* — Vc~ SV* 
that verifies the condition 



disv*) 



(9) 



From the above considerations, we conclude that cooling 
of the nanotube vibrations can only occur for bias volt- 
ages below the Coulomb blockade threshold [SV > 0). 

However, below the Coulomb blockade threshold volt- 
age, charge transport is blocked at zero temperature. The 
temperature needed to overcome the Coulomb blockade is 
determined by fcsT > eSV . On the other hand, the tem- 
perature cannot be too high, since otherwise "backward" 
transitions from the leads to the nanotube would no 
longer be negligible and possibly compensate for the vi- 
brons absorbed during the "forward" transitions. These 
conditions restrict the range of possible temperatures to 
the interval: eVc > ksT > eSV eSV* . The or- 
der of magnitude of the lower bound can be found by 
means of the condition ([9]) and by estimating the ca- 
pacitance between the STM tip and the nanotube as 
Cs — 27TeQD / \n{2h/ro), where h ^ 1, Eq is the vacuum 
permittivity, I? ^ 10 nm is the characteristic diameter of 
the STM tip, tq ^ 0.5 nm is the nanotube radius. One 
finds that the temperature required in order to overcome 
the Coulomb blockade at SV* is about 0.1 K. 

For a quantitative analysis of the cooling mechanism 
described above, we followed the standard procedure to 
derive a generalized master equation for the reduced den- 
sity matrix that describes the nanotube degrees of free- 
dom . After tracing out the charge degrees of freedom 
and applying a perturbation approach with respect to the 
small parameters et.p one gets a set of equations for the 
probabilities p„ to find the nanotube in the Fock state 
\n) characterized by n vibrons. If the rate F^ of tunnel- 
ing from the nanotube to the leads is much larger than 
the rate Ts of tunneling from the STM to the nanotube 
these equations reduce to 



= (10) 



[{An + 2)4 + {2n + 1)2^ _ 2spet] Pn 
[{Sp + £tf + £p] (n -I- l)pn+l + [(£p - etf + £p] npn-i , 
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Figure 2. The average number (n) of vibrons is plotted against 
the difference = Vc — V between the Coulomb blockade 
threshold voltage Vc and the bias voltage V . Each curve cor- 
responds to a different quality factor of the oscillator, while 
the straight line gives the thermal average number of vibrons 
at the temperature of T = 1 K. 

where C-y describes the interaction with the environ- 
ment, which takes the standard form [l^: Cj[pn] = 

j{n+l)[{nth + l)pn+i~nthPn]-7'n[{nth + i)Pn-nthPn-i], 
where 7 = cu/Q, Q being the quality factor of the nan- 
otube resonator, and nth = (e''"/'^^^— 1)"-'^ is the thermal 
average number of vibrons. 

Equation pH)) can be solved for the stationary proba- 
bility distribution p„ with the result 

p„ = (1 - r)r" , (11) 
_ e'p + jet - £pf + (7/r5)"-th 

)' + (7/r5)Kh + i) 

The average number of vibrons, (n) — ^^^^ mpm — 
r/(l — r) is plotted as a function of the bias voltage 
for several values of the quality factor in Fig. [51 The 
theoretical limit of the cooling efficiency occurs for the 
bias voltage defined by ep[5V*) = et/^/2) in the limit 
Q ^ 00. Equation (jlll) implies that the corresponding 
average number of excitations is nmin = (v^— 1)/2 ~0.2. 

In order to investigate the signatures of the cooling 
mechanism in a directly measurable property, we have 
calculated the current / perturbatively to second order 
in Et^p with the result 

/ = /o(l+£?(l + 2(n))) . (12) 

Here /q = er5(V^,T) with Tg ~ hgT if hgT > e5V* 
and Ts remains independent of voltage in a certain 
voltage interval, where the differential conductance will 
be completely determined by the derivative of the av- 
erage number of vibrons with respect to voltage, i.e. 



81/ dV = 2/oe^9(n)/9V. Therefore, the cooling effect 
will be reflected in the structure of the dl/dV — V curves 
and accessible for experimental investigation. 

In conclusion we have proposed a novel mechanism for 
ground-state cooling of nanomechanical resonators based 
on the injection of a tunneling current from a voltage- 
biased STM tip. For the model system considered we 
have shown that the direction of the electrostatic force 
that acts on a suspended charged nanotube can be flipped 
by changing the voltage bias. This makes it possible to 
control the interference between two distinct contribu- 
tions to the quantum mechanical probability amplitudes 
for vibron absorption and emission during electron tun- 
neling. For a bias voltage slightly below the Coulomb 
blockade threshold voltage, the probability amplitude for 
vibron emission becomes very small. At this bias a ther- 
mally activated current therefore leads to a cooling of 
the nanomechanical vibrations. Our analysis shows that 
the effective temperature that can be reached may corre- 
spond to an average vibron population of the fundamen- 
tal bending mode as low as 0.2. The cooling mechanism, 
which should be observable by its effect on the differ- 
ential conductance of the system, is crucially dependent 
on the Coulomb blockade phenomenon and hence on the 
quantization of electric charge. 

Partial support from the Swedish VR and SSF, the 
EC project QNEMS (FP7-ICT-233952), the Faculty of 
Science at the University of Gothenburg through its 
Nanoparticle Research Platform and from the Korean 
WCU program funded by MEST/NFR (R3 1-2008-000- 
10057-0) is gratefully acknowledged. 



[1] A. D. O'Connell et al.. Nature 464, 697 (2010). 

[2] M. Blencowe, Phys. Rep. 395, 159 (2004). 

[3] 1. Wilson-Rae, P. ZoUer, and A. Imamoglu, Phys. Rev. 
Lett. 92, 075507 (2004). 

[4] S. ZippiUi, G. Morigi, and A. Bachtold, Phys. Rev. Lett. 
102, 096804 (2009); F. Pistolesi, J. Low Temp. Phys. 
154, 199 (2009); S.-H. Ouyang, J. Q. You, and F. Nori, 
Phys. Rev. B 79, 075304 (2009); G. Sonne et al., Phys. 
Rev. Lett. 104, 226802 (2010). 

[5] B. J. LeRoy, S. G. Lemay, J. Kong, and C. Dekker, Na- 
ture 432, 371 (2004). 

[6] L. D. Landau and E. M. Lifshitz, Electrodynamics of con- 
tinuous media (Butterworth-Heinemann, Oxford, 1998). 

[7] L. M. Jonsson, L. Y. Gorelik, R. I. Shekhter, and M. Jon- 
son, Nano Lett. 5, 1165 (2005). 

[8] T. Nord and L. Y. Gorelik, Low Temp. Phys. 31, 534 
(2005). 

[9] L. Y. Gorelik, D. Fedorets, R. I. Shekhter, and M. Jonson, 
New J. Phys. 7, 242 (2005). 
[10] H.-P. Breuer and F. Petruccione, The theory of open 
quantum systems (Oxford University Press, 2002). 



